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ABSTRACT 

Conventionally, long GRBs are thought to be caused by the core collapses of massive 
stars. During the lifetime of a massive star, a stellar wind bubble environment should 
be produced. Furthermore, the microphysics shock parameters may vary along with 
the evolution of the fireball. Here we investigate the variation of the microphysics 
shock parameters under the condition of wind bubble environment, and allow the 
microphysics shock parameters to be discontinuous at shocks in the am- 
bient medium. It is found that our model can acceptably reproduce the 
rebrightenings observed in GRB afterglows, at least in some cases. The 
effects of various model parameters on rebrightenings are investigated. The rebright- 
enings observed in both the R-band and X-ray afterglow light curves of GRB 060206, 
GRB 070311 and GRB 071010 A are reproduced in this model. 
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1 INTRODUCTION 

Gamma-ray bursts (GRBs) are attractive astrophysics phe- 
nomena and had puzzled astronomers for about twenty- 
four years after their discovery in 1973 (Klebesadel et al. 
1973). The discovery of long-lived, multi-band counterparts 
of GRBs, namely, the afterglows of GRBs, in 1997 is a wa- 
tershed in GRB research (Costa et al. 1997; Van Paradijs et 
al. 1997; Frail et al. 1997). Soon after that, the so-called fire- 
ball model is recognized as the standard model in view of the 
fact that it can explain most features of GRB observations 
well. However, as the advance of observation techniques and 
the accumulation of observational data, especially after the 
launch of Swift satellite (Gehrels 2004), a lot of unexpected 
behaviours appear in GRB afterglows, such as the canonical 
steep- shallow- normal decay and flares in X-ray afterglows, 
the flattish decay phase and various rebrightenings in optical 
afterglows (For review, see Zhang 2007). 

In fact, there are more and more rebrightenings detected 
in GRB afterglows, including GRB 970508 (Galama et al. 
1998a), GRB 990123 (Sari & Piran 1999), GRB 021004 (Laz- 
zati et al. 2002), GRB 030329 (Berger et al. 2003), GRB 
050525A (Blustin et al. 2006), GRB 050820A (Cenko et al. 
2006), GRB 050721 (Antonelli et al. 2006), GRB 060206 
(Stanek et al. 2007; Liu et al. 2008), GRB 070125 (Up- 
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dike et al. 2008), GRB 070311 (Guidorzi et al. 2007a), GRB 
071003 (Perley et al. 2008), GRB 071010A (Covino et al. 
2008). Many different mechanisms have been proposed to 
explain these rebrightenings, such as density jump (Lazzati 
et al. 2002; Dai k Wu 2003; Tam et al. 2005), energy in- 
jection (Huang et al. 2006), two-component jet (Huang et 
al. 2004, 2006; Liu et al. 2008), reverse shock (Sari & Pi- 
ran 1999), reverberation of the energy input measured by 
prompt emission (Vestrand et al. 2006), turn-on of the ex- 
ternal shock (Stanek et al. 2007; Molinari et al. 2007), large 
angle emission (Panaitescu & Kumar 2007), and spectral 
peak of existing forward shock (Shao & Dai 2005). 

Among all the mechanisms, the density jump model 
needs to be paid special attention to, because density jump 
in surrounding medium of GRBs is a very natural hy- 
pothesis. Since there are more and more examples indicat- 
ing that some GRBs are associated with Type Ic super- 
novae (e.g. SN 1998bw/GRB980425, Galama et al. 1998b; 
SN2003dh/GRB030329, Price et al. 2003) and many host 
galaxies of GRBs are in process of active star formation 
(Fruchter et al. 1999; Djorgovski et al. 1998), it is be- 
lieved that the progenitors of long GRBs are massive Wolf- 
Rayet (WR) stars (Woosley 1993). Massive stars usually 
produce very strong stellar wind to push the initial inter- 
stellar medium (ISM) in their neighborhood away during 
their lifetime. The surrounding of these GRBs should be a 
wind bubble following the density profile p oc rather 
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than the usual homogeneous ISM. Several authors (Castor 
et al. 1975; Weaver 1977; Ramirez-Riuz et al. 2001) further 
found that beyond some typical radius of the wind bubble, 
the swept- up mass is too large to be pushed by the wind. As 
a result, the wind materials pile up at the edge of the wind 
bubble to form a density jump. 

Lazzati et al. (2002) proposed that a density jump can 
lead to a rebrightening in the afterglow. In usual case, when 
the observing frequence is between the peak frequence (z^m) 
and the cooling frequence (z^c), the amplitude of the rebright- 
ening should be proportional to the square root of the den- 
sity contrast. However, if the density contrast is too high, 
the rebrightening will be weakened, since Vc will decrease 
and become less than the observing frequence. In the more 
detailed numerical simulations by Huang et al. (2007), no ob- 
vious rebrightenings can be discriminated when the density 
jump amplitude is set to 100, which suggests that the simple 
density jump model is not an ideal mechanism to produce 
the rebrightenings. In the recent studies by Nakar & Granot 
(2007) and van Eerten et al. (2009), a full treatment of the 
transient features at the jump moment, even including the 
thickness of the blastwave, was considered. Thus these stud- 
ies should be a more accurate approximation to the reality. 
Interestingly, no obvious rebrightenings associated with the 
density jumps are found. 

In short, a simple density jump model is difficult to 
explain the observed rebrightenings in GRB afterglows. On 
the other hand, the early afterglows of most GRBs exhibit 
flattish decays with a < 0.8, where a is defined as Fi, oc 
t~^. This is difficult to explain using the standard fireball 
model when the density profile is p oc . Some previous 
works have suggested that the microphysics parameters may 
vary during the evolution of the ffieball (Rossi & Rees 2003; 
loka et al. 2006; Fan & Piran 2006; Panaitescu et al. 2006; 
Granot, Konigl & Piran 2006). We believe that this kind of 
variation could resolve these problems. 

In this paper, we show that the observed rebrighten- 
ings in GRB afterglows can be well reproduced by assuming 
varying microphysics shock parameters associated with the 
wind bubble environments. The outline of our paper is as 
follows: in §2 we introduce our model detailedly. We then 
numerically investigate the effects of various parameters on 
the afterglows in §3, and reproduce the R-band and X-ray af- 
terglow light curves of GRB 060206, GRB 070311 and GRB 
071010A in §4. Our discussion and conclusions are presented 
in §5. We use an assumptive cosmology of Hq = 65 km s~^ 
Mpc~^, Qm = 0.30 and Qa = 0.70 throughout the paper. 



2 MODEL 

2.1 Dynamics and radiation process 

In the standard fireball model, the outflow of GRB, 
which moves relativistically, interacts with the surrounding 
medium to form an external shock. A constant fraction 
of the shock energy will be transferred to the swept-up elec- 
trons and accelerate them to relativistic velocity. Similarly, a 
constant fraction of the shock energy will go to the mag- 
netic field. These shocked relativistic electrons move in the 
magnetic field and emit synchrotron radiation to produce 
broadband afterglows. 



We use the convenient equations developed by Huang 
et al. (1999, 2000a, 2000b, 2003) to describe the dynamics 
and radiation process of the eject a. The evolution of the 
bulk Lorentz factor 7, the shock radius and the swept-up 
medium mass m, is described by three differential equations. 



: /3c7(7+ vV-1), 



(1) 



(2) 



(3) 



d7 7 ~ 1 

dm Mej + em + 2(1 — e)jm ' 

dR 
dt 

dR 

where nip is proton mass, Mej is the initial mass of the 
outflow, is the half opening angle of the jet, n is the 
number density of the environment, /3 — y^t^ — I/7, and 
7 ^ (47 + 1)7(37) is the adiabatic index (Dai, Huang, & 
Lu 1999). e is the radiative efficiency, which equals 1 for a 
highly radiative case, and equals in the adiabatic case. We 
ignore the sideways expansion of the jet in our model, be- 
cause many numerical simulations indicate that it is a very 
slow process (Granot et al. 2001; Cannizzo, Gehrels & Vish- 
niac 2004; Zhang & MacFadyen 2009). We consider that the 
afterglow flux mainly comes from the synchrotron radiation 
of the shocked relativistic electrons. 

Since the speed of light ( c ) is finite, photons received by 
the observer at a particular time t are not radiated simul- 
taneously, but come from a distorted ellipsoid determined 
by 



t ■■ 



r- 



pcosO 



dR = const, 



(4) 



within the jet boundaries, where G is the angle between the 
velocity of emitting material and the line of sight. This is the 
so-called EATS effect. We integrate over the EATS to obtain 
an accurate observed flux in our simulations (see Huang et 
al. 2007 for more details). 

2.2 Environment 

Massive stars usually produce very strong stellar wind dur- 
ing their lifetime. This stellar wind interacts with the initial 
ISM and forms two shocks: a reverse shock propagates back 
into the stellar wind and a forward shock that propagates 
into the ISM. The resulting surrounding medium is broken 
into four regions (Castor et al. 1975; Weaver 1977; Ramirez- 
Riuz et al. 2005; Pe'er and Wijers 2006) as shown in Figure 1. 
They are (from inside to outside): (1) the unshocked stellar 
wind; (2) the shocked stellar wind; (3) the shocked ISM; (4) 
the unshocked ISM. 

Region (1) is a stellar wind bubble. The number density 
of this region is given by 



ni 



M 



(5) 



where M is the mass loss rate of the progenitor, t;w is the 
launching speed of the stellar wind and r is the radial dis- 
tance from the central star. We fix v^^ as 10^ cm s~^ and 
only vary M in our simulations, since the density profile of 
the environment only depends on the ratio M /v^. 

As the stellar wind expands, the swept-up mass of the 
ISM increases and is finally comparable to the mass in the 
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Radius from progenitor 

Figure 1. The illustration of the density profile used in our 
model. 



wind at a radius rwind, which is determined by (Ramirez- 
Ruiz et al. 2001, 2005) 



"^wind 



1.6 X io'^'m; 



18 A>3/10^ 



-3/10 1/10,2/5 
c:i\/r '^w s cm, 



(6) 



where nisM,3 is the initial density of the homogeneous ISM, 
in units of 10^ cm~^, and t^,6 is the life time of the WR 
phase of the star, in units of 10^ yr. Beyond this radius, 
in region (2), the swept-up mass is too large to be pushed 
away by the wind. The material piles up in this region to 
form a density jump. Because of the effect of the reverse 
wind shock, the material is hot in this region. The density 
is approximately constant and equals to n2 = 4ni(rwind), 
where ni(rwind) is the density in region (1) at the radius of 

^ — "^wind- 

As shown by Castor et al. (1975), the boundary between 
Region (2) and Region (3) is at the radius of 



r2-3 



1.6 X 10"My/n-;/>:(|t!f|cm. 



(7) 



Region (3) is a thin, dense shell containing most of the 
swept-up ISM. Its width is only about r2-3/12 and the den- 
sity in this region is about 4nisM- Here nisM is also the 
density in region (4). 

We use the density profile introduced above as the en- 
vironment surrounding GRBs in this paper. As discussed by 
Pe'er and Wijers (2006), for a GRB event, the blast wave can 
not reach Region (3) in the relativistic phase. Our numeri- 
cal results also prove that the fireball is in Region (1) and 
Region (2) during typical observable time. So we only use 
Region (1) and Region (2) as the environment in our work. 
For simplicity, we take rwind as a free parameter in our cal- 
culations, so that we do not need to consider the detailed 
values of and t^. 

Unfortunately, the simple density jump model alone is 
difficult to explain the rebrightenings in GRB afterglows 
(Huang et al. 2006, 2007; Nakar & Granot 2007). Moreover, 
the early optical and/or X-ray afterglow light curves of some 
GRBs are flatter than the prediction of a normal fireball in 
a simple wind environment. So we need to consider other 
effects. Varying microphysics shock parameters should be a 
possible solution. 



2.3 Microphysics parameters 

In the standard afterglow model, it is usually assumed that 
the microphysical parameters are constant throughout the 
evolution of the fireball. However, the actual microphysical 
processes in the relativistic shocks, for example, the energy 
transfer from protons to electrons and magnetic fields, are 
still not clearly known. It is possible that the microphysi- 
cal parameters are variational. In fact, the variation of the 
microphysics parameters during the evolution of the fireball 
has been considered in many previous studies (Rossi & Rees 
2003; loka et al. 2006; Fan & Piran 2006; Panaitescu et al. 
2006; Granot, Konigl k Piran 2006). Fan & Piran (2006) 
and Panaitescu et al. (2006) have engaged the assumption 
that the equipartition factors and are functions of 7 to 
explain the shallow decay phase in some X-ray afterglows. 
Here we will use a similar idea to investigate the observed 
rebrightenings in GRB afterglows. 

In our model, the circum-burst environment is divided 
into two parts, i.e. Region (1) and Region (2). The material 
in Region (1) is the unshocked stellar wind thrown by the 
mass star. In Region (2) the material consists of the shocked 
stellar wind mixed with a small fraction of the swept-up 
ISM. It is hot because of the effect of the reverse wind shock 
(Castor et al. 1975; Weaver 1977; Ramirez-Riuz et al. 2005; 
Pe'er and Wijers 2006). We can imagine that the physical 
condition, such as the strength of the magnetic field, the 
temperature and density of the material, could be different 
between these two regions, so the evolution of microphysics 
parameters may also be different accordingly. In our study, 
we use different parameters for these two regions. We assume 
that 

ee = 6,o7""\ (8) 

= ^5,07""^^ (9) 

in Region (1), and 

= 6,07""', (10) 
= ^5,07""^' (11) 

in Region (2). 



3 NUMERICAL RESULTS 

In this section, we present our numerical results to show the 
effects of various physical parameters on the R-band light 
curves, using the model described in §2. For convenience, 
we first define a set of "standard" parameters, as shown in 
the second column of Table 1. 

Figure 2 illustrates the effect of the parameter Oo on the 
R-band light curves. The solid line corresponds to the "stan- 
dard" parameters. The dashed line corresponds to Oo — 0.3 
rad and the dotted line corresponds to Oo — 0.03 rad. We 
can see that when the change of microphysics parameters is 
considered, an obvious rebrightening appears. Interestingly, 
when Oo is larger, the duration of the rebrightening becomes 
longer. This is not difficult to understand. The afterglow 
brightness is dominated by the high latitude emission, so 
that the EATS shows a ring- like pattern. As a result, the 
brightness will be kept on a relatively higher level after the 
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Table 1. Modeling Parameters 



Parameters 


Standard 


GRB 060206 


GRB 070311 


GRB 071010A 


z 


1.0 


4.045 


1.0 


0.98 


Eo (ergs) 


1.0 X 10^3 


9.0 X 10^2 


8.0 X 10^1 


1.3 X 10^2 


Oo (rad) 


0.1 


0.06 


0.08 


0.11 


M (Mo/yr) 


1.0 X 10-6 


2.0 X 10-^ 


6.0 X 10-6 


4.0 X 10-^ 


^wind (cm) 


1.0 X lO^s 


3.0 X 10^6 


2.0 X 10^'^ 


5.1 X 10^6 


p 


2.5 


2.2 


2.3 


2.1 




0.1 


0.3 


0.3 


0.3 




0.01 


0.03 


0.03 


0.03 


al 


1.5 


1.2 


0.9 


1.5 


a2 


1.3 


0.6 


0.0 


0.6 




10' 10' 10" 

Time Since Burst (s) 




10' 10' 10" 

Time Since Burst (s) 



Figure 2. The effect of the parameter on the R-band hght 
curve. The sohd hne corresponds to the standard parameters. 
The dashed hne corresponds to Oq = 0.3 rad and the dotted hne 
corresponds to 6o = 0.03 rad. 



Figure 3. The effect of the parameter M on the R-band hght 
curve. The sohd hne corresponds to the "standard" parameters. 
The dashed hne corresponds to M = 1.0 x 10—^ Mq /yr and the 
dotted hne corresponds to M = 1.0 X 10-^ Mq/yt. 



density jump moment, to form a flat phase. When Oq be- 
comes larger, the time delay of the photons from the edge of 
the outflow also becomes larger. So the flat phase is longer. 

Environment is an important factor that affects the evo- 
lution of GRB afterglows. As mentioned before, the charac- 
ter of wind environment is mainly determined by the ratio 
M/vw We fix as 10^ cm s~^ and only vary M here to 
investigate the effect of the parameter M on the R-band 
light curves. The results are shown in Figure 3. The solid 
line corresponds to the "standard" parameters. The dashed 
line corresponds to M = 1.0 x 10~^ Mq/jt and the dotted 
line corresponds to M = 1.0 x 10~^ Mq/yt. It is clearly 
seen that for a smaller M value, the rebrightening appears 
earlier and the flux before the rebrightening are lower. This 
result is easy to understand. A smaller M corresponds to a 
smaller density, and the deceleration of the external shock 
is slower. So the Lorentz factor is larger and the fireball can 
meet the density jump earlier. A larger Lorentz factor also 
makes and smaller and decreases the flux. 

The parameter rwind is another factor to describe the 
character of the environment. The effect of rwind on the R- 
band light curve is illustrated in Figure 4. The solid line 
corresponds to the "standard" parameters. The dashed line 



corresponds to rwind = 5.0 x lO^^cm and the dotted line 
corresponds to rwind = 5.0 x lO^^cm. We can see that the 
rebrightening appears earlier when rwind is smaller, because 
rwind determines the position of the density jump. 

The effect of the parameter al on the R-band light 
curve is shown in Figure 5. The solid line corresponds to 
the "standard" parameters. The dashed line corresponds to 
al = 1.7 and the dotted line corresponds to al = 1.3. It is 
obvious that a larger al makes the flux before the rebright- 
ening lower. This is not difficult to understand. With the 
increase of al, the values of and before the rebright- 
ening become smaller. This suppresses the radiation flux at 
early time. 

The effect of the parameter a2 on the R-band light 
curve is illustrated in Figure 6. The solid line corresponds 
to the "standard" parameters. The dashed line corresponds 
to a2 = 1.5 and the dotted line corresponds to a2 = 1.1. 
As expected, the effect of a2 is similar to that of al for the 
same reason. 
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Figure 4. The effect of the parameter rwind on the R-band hght 
curve. The sohd hne corresponds to the "standard" parameters. 
The dashed hne corresponds to rwind = 5.0 x lO-'^^cm and the 
dotted hne corresponds to rwind = 50 x lO-'^^cm. 
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Figure 5. The effect of the parameter al on the R-band hght 
curve. The sohd hne corresponds to the "standard" parameters. 
The dashed hne corresponds to al = 1.7 and the dotted hne 
corresponds to al = 1.3. 



4 COMPARISON WITH OBSERVATIONS 

In this section we select three GRBs with obvious rebright- 
enings in both optical and X-ray bands, i.e. GRB 060206, 
GRB 070311 and GRB 071010A, and reproduce their R- 
band and X-ray afterglow light curves numerically, using 
the model described in §2. 



4.1 GRB 060206 

GRB 060206 was detected at 04:46:53 UT on Feb. 6 by 
the Burst Alert Telescope (BAT) onboard the Swift satel- 
lite (Morris et al. 2006). The temporal profile of the BAT 
light curve shows a single peak, with a total fluence of 
8.4 ± 0.4 X 10"'^ erg cm"^ in the 15-350 keV band (Mor- 
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Figure 6. The effect of the parameter a2 on the R-band light 
curve. The solid line corresponds to the "standard" parameters. 
The dashed line corresponds to a2 = 1.5. and the dotted line 
corresponds to a2 = 1.1. 



ris et al. 2006; Palmer et al. 2006). The total duration of the 
burst is Tgo = 7 ± 2s, suggesting that it is a long-duration 
GRB (Palmer et al. 2006). 

The UVOT telescope onboard the Swift satellite first 
observed the optical afterglow of GRB 060206 and found 
it located at the position of a = 13^31'^43^416, 6 = 
+35°03'03''.6 (J2000.0; Boyd et al. 2006), corresponding 
to the Galactic extinction of Ar = 0.033 mag (Schlegel, 
Finkbeiner & Davis 1998). After that, a number of tele- 
scopes made detailed follow up observations, such as the 2- 
m robotic Liverpool Telescope (Monfardini et al. 2006), the 
Rapid Telescopes for Optical Response (RAPTOR; Wozniak 
et al. 2006) and the MDM telescope (Stanek et al. 2007). 
From the calibrated ALFOSC spectrum of the afterglow, 
Fynbo et al. (2006) determined the redshift of GRB 060206 
as z = 4.045. X-ray afterglow was also detected by XRT on 
Swift at 58 s after the trigger (Morris et al. 2006) . 

There are significant rebrightenings in both the optical 
and X-ray afterglows of GRB 060206 at about 3000 s after 
the burst. Liu et al. (2008) suggested that the rebrightening 
comes from a jet-like ejecta with a larger viewing angle, 
which is produced by the central engine. 

Here we try to use our model described in §2 to re- 
produce the rebrightenings in both R-band and X-ray light 
curves of GRB 060206. Our best-fit physical parameters are 
presented in Table 1, and our modeling curves are illustrated 
in Figure 7. The observed R-band data are taken from Mon- 
fardini et al. (2006), Wozniak et al. 2006 and Stanek et al. 
(2007), and the X-ray data are taken from the Swift XRT 
light curve repository^ (Evans et al. 2007). We can see in 
Figure 7 that our model can generally reproduce the light 
curves well. Note that at the late stage, the observed X-ray 
light curve is obviously too flat as compared with our theo- 
retical result. This may be due to the contamination from a 
nearby source, as already mentioned by Stanek et al. (2007). 
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Figure 7. Our best fit to the R-band and X-ray afterglow light 
curves of GRB 060206. The lines are our theoretical curves and 
the fitting parameter values are given in Table 1. The R-band 
observed data are taken from Monfardini et al. (2006), (Wozniak 
et al. 2006) and Stanek et al. (2007), and the X-ray data are taken 
from the Swift XRT light curve repository (Evans et al. 2007). All 
optical data points have been corrected for the Galactic extinction 
(Schlegel, Finkbeiner & Davis 1998). 
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Figure 8. Our best fit to the late rebrightenings in the R-band 
and X-ray light curves of GRB 070311. The lines are our theoret- 
ical curves and the fitting parameter values are given in Table 1. 
The R-band observed data are taken from Guidorzi et al. (2007a; 
for REM data) and GCN circulars (Cenko 2007; Dai et al. 2007; 
Garnavich et al. 2007; Halpern &: Armstrong 2007a, b,c,d; Jelmek, 
Kubanek & Prouza 2007; Kann 2007; Wren et al. 2007), and the 
X-ray data are taken from the Swift XRT light curve repository 
(Evans et al. 2007). All optical data points have been corrected 
for the Galactic extinction (Schlegel, Finkbeiner & Davis 1998). 



4.2 GRB 070311 

GRB 070311 occurred at 01:52:50 UT on 2007 March 11 and 
was detected by the INTEGRAL/IBAS in IBIS/ISGRI data. 
The duration of the 7-ray prompt emission is about 50 s in 
the 20-200 keV band, with the peak flux of 0.9 photons cm~^ 
s""*^ (1 s integration time) and the total fluence of 2 x 10~^ 
erg cm~^ (Mereghetti et al. 2007). 

The afterglow of GRB 070311 was first discovered 
by the Rapid Eye Mount (REM) telescope at 55 s af- 
ter the burst. A bright fading source was found at a = 
05^50^08^.21, 5 = +03°22'30''.3 (J2000.0; Covino et al. 
2007). This location has a low Galactic latitude (/ = 
202°. 766, 6 = -11°.998 in Galactic coordinates), so the 
Galactic extinction is as much as Ar — 2.038 mag (Schlegel, 
Finkbeiner & Davis 1998). X-ray afterglow was also found 
by the observations of XRT on Swift at about 7000 s after 
the trigger (Guidorzi et al., 2007b). 

The optical afterglow of GRB 070311 has complex 
structures. Besides a normal power law decay, there are two 
fast rise exponential decay (FRED) shaped pulses peaking 
around 80 and 140 s after the peak of the GRB and possibly 
accompanied by the tail of prompt 7-ray emission (Guidorzi 
et al. 2007a). Another structure, which is the most attractive 
feature, is a significant late rebrightening between 3 x 10^ 
s and 2 X 10^ s after the trigger in both optical and X-ray 
bands (Guidorzi et al. 2007a). Guidorzi et al. (2007a) have 
used a power law function plus a FRED shaped pulse to 
fit the afterglows in both R and X-ray bands. More phys- 
ically, they suggested that the late afterglow rebrightening 
of GRB 070311 can come from a refreshed shock. Addition- 
ally, they argued that the density jump in the surrounding 
medium can be the origin of the rebrightening too, but less 
appealing. 

Here we use our model described in §2 to reproduce the 
late rebrightenings in both R-band and X-ray light curves 



of GRB 070311. We take the observed R-band data from 
Guidorzi et al. (2007a; for REM data) and GCN circulars 
(Cenko 2007; Dai et al. 2007; Garnavich et al. 2007; Halpern 
& Armstrong 2007a,b,c,d; Jelmek, Kubanek & Prouza 2007; 
Kann 2007; Wren et al. 2007), and the X-ray data from 
the Swift XRT light curve repository (Evans et al. 2007). 
Because the redshift of GRB 070311 is still unknown, we 
simply assume a redshift of 2; = 1.0 in our calculation. Our 
best-fit physical parameters are presented in Table 1, and 
our modeling curves are illustrated in Figure 8. We see that 
the interpretation of our model to this event is also generally 
acceptable. 

4.3 GRB 071010A 

GRB 071010A triggered the Swift satellite at 03:41:12 UT 
on 2007 Oct. 10 (Moretti et al. 2007). It is a traditional long 
GRB with the burst duration of T90 = 6 d= Is. The total 
fluence of the burst is about 2.0 x 10~^ erg cm~^ (Krimm 
et al. 2007). 

TAROT telescope identified the optical afterglow of 
GRB 071010A firstly, only 124 s after the trigger (Klotz, 
Boer & Atteia, 2007). Bloom et al. (2007) further found 
that the position of the burst was a ^ 19^12'^14\624, 
6 = -32°24'0rM6 (J2000.0), corresponding to the Galac- 
tic extinction of Ar — 0.263 mag (Schlegel, Finkbeiner & 
Davis 1998). Spectral observations were done subsequently 
by Keck and the analysis of the Mg H and Fe H lines gave 
the redshift of the burst as z = 0.98 (Prochaska et al. 2007). 
X-ray data were also acquired by XRT observations on Swift 
(Guidorzi et al. 2007b), but no radio afterglow was detected 
at 8.46 GHz until almost 2 days after the burst (Chandra & 
Frail 2007). 

The most interesting feature of the afterglows of GRB 
071010A is the appearance of a sharp rebrightening at about 
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Figure 9. Our best fit to the R-band and X-ray afterglow light 
curves of GRB 071010A. The lines are our theoretical curves and 
the fitting parameter values are given in Table 1. The R-band 
observed data are taken from Covino et al. (2008), and the X-ray 
data are taken from the Swift XRT light curve repository (Evans 
et al. 2007). All data points have been corrected for the Galactic 
extinction (Schlegel, Finkbeiner & Davis 1998). 



0.6 days after the GRB, clearly seen in both optical and 
X-ray bands. Covino et al. (2008) have used a Beuermann 
et al. (1999) function, and a simple step function for the 
rebrightening, to fit the afterglows. They suggested that the 
rebrightening was due to an injection of abundant energy, 
which was comparable to the initial energy in the fireball. 
The steepening after the rebrightening was interpreted as a 
jet break at around 1 day after the burst. 

We try to use our model introduced in §2 to reproduce 
the R-band and X-ray light curves of GRB 071010A. In our 
modeling, the observed R-band data are taken from Covino 
et al. (2008), and the X-ray data are taken from the Swift 
XRT light curve repository (Evans et al. 2007). Our best- 
fitting physical parameters are presented in Table 1, and our 
modeling curves are illustrated in Figure 9. It is clear that 
our model can reproduce the rebrightening features of GRB 
071010A very well. 



5 DISCUSSION AND CONCLUSION 

In the standard fireball model, it is assumed that the circum- 
burst medium density is constant or following a single p oc 

law. The microphysics shock parameters are also usu- 
ally assumed invariable during the evolution of the fireball. 
These models can basically explaining many pre- Swift GRB 
afterglows, whose spectrum and light curves can be approxi- 
mated as broken power-law functions (Panaitescu & Kumar 
2001a, 2001b; Yost et al. 2003). The launch of Swift satel- 
lite (Gehrels 2004) makes it possible to observe early after- 
glows of many GRBs in the first few hours after the trigger. 
Many remarkable and unexpected features are found, such 
as marked rebrightenings and the flattish decay phase in 
early GRB afterglows. 

In fact, the standard fireball model is obviously too sim- 
plified. If long GRBs indeed originate from the death of mas- 
sive stars, the environment of GRBs should have complex 



structures rather than have a constant or single power law 
profile. Some analytical (Castor et al. 1975; Weaver 1977; 
Ramirez-Riuz et al. 2005; Pe'er and Wijers 2006) and nu- 
merical (Ramirez- Riuz et al. 2001) studies suggest that the 
circum-burst density profile should be a wind bubble, as- 
sociated with a few density jumps. On the other hand, the 
microphysics parameters, such as and ^b, may vary dur- 
ing the evolution of the fireball (Rossi & Rees 2003; loka 
et al. 2006; Fan & Piran 2006; Panaitescu et al. 2006; Gra- 
not, Konigl & Piran 2006). Actually, the fast decrease of the 
cooling frequency Uc in some GRBs suggests that may 
be evolving (Panaitescu et al. 2006). Some previous studies 
also suggest that and ^b may be different for the forward 
shock and the reverse shock (Fan et al. 2002; Wei, Yan & 
Fan 2006), and may also be different for Region (1) and Re- 
gion (2) (Gebdre et al. 2007; Kamble, Resmi & Misra 2007). 
So, and ^b may depend on the strength of the shock and 
the environment. 

In our study, we combine the wind bubble environments 
and the change of microphysics shock parameters together. 
Comparing with the standard fireball model, our model has 
three more parameters (i.e. Twind, Oil and a2). We find that 
this model can produce the observed rebrightenings and flat- 
tish decay in GRB afterglows successfully. We have shown 
that the observed rebrightenings in both the optical and X- 
ray afterglow light curves of GRB 060206, GRB 070311 and 
GRB 071010A can be well explained by our model. 

We have also investigated the effects of various param- 
eters on the light curves numerically. We can imagine that if 
the stellar wind produced by the progenitor is very strong, 
or the launching speed of the stellar wind t'w is small, or 
the value of rwind is large enough, there will be no rebright- 
ening within the usual observation time. So, basically the 
wind bubble environment can also give birth to a steadily 
decaying afterglow that shows no rebrightening. 

In our work, we neglect the effect of the reverse shock. 
According to the study by Dai & Lu (2002), when an ultra- 
relativistic blast wave interacts with a density jump medium, 
the resulting reverse shock is relativistic only if the ampli- 
tude of the density jump is much larger than 21. In our 
model, the amplitude of the density jump is only 4, so 
the corresponding reverse shock is Newtonian. The emis- 
sion from the Newtonian reverse shock is very weak, and 
can be omitted. Moreover, although we use an abrupt den- 
sity jump, the actual increase of density may be gradual. In 
this situation, the emission from the reverse shock will even 
be much weaker. 

Currently, a complete understanding of the microphys- 
ical processes in the relativistic shocks is still lacking. Using 
the derived microphysical parameters from GRB modeling, 
people may be able to get some constraints on the shock 
physics. The derived parameter values of M and nisM are 
also very important. They are closely related to the evo- 
lution and the environment of the massive star. They may 
give some hints on the characters of the progenitor, such as 
the initial main-sequence mass and the metallicity. So, from 
the modeling parameters, we can also obtain some useful 
information about GRB origin. 
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